ON COMPARING THE WRITHE OF A SMOOTH CURVE 
c ^". TO THE WRITHE OF AN INSCRIBED POLYGON 

O 

^ ! JASON CANTARELLA 

D ' Abstract. We find bounds on the difference between the writhing numbers of a smooth curve 

and a polygonal curve inscribed within. The proof is based on an extension of Fuller's difference 
£NJ ' of writhe formula to the case of polygonal curves. The results establish error bounds useful in the 

, numerical computation of writhe. 



1. Introduction 



o 

Q 

-s, 

The writhing number measures the wrapping and coiling of space curves. Writhe has proved 
useful in molecular biology, where it is used to study the geometry of tangled strands of DNA 
[|T7p; often with the famous Calugareanu-White formula for a curve C in space with a normal 
field V [00^0]: 

Lk{C,C + eV) = Tw(C,V) +Wr(C). 

m ■ 

In these applications, and in numerical simulations performed by biologists and mathematicians, it 
CnI \ is often required to compute writhing numbers using numerical methods. 

^ | Several authors have presented algorithms for computing the exact writhing number of an pi- 

edge polygonal curve in a finite number of steps [H |], [T], \V7\\ . The fastest of these algorithms runs 
r-| ■ in time 0(n L6 ), while earlier methods use time 0(n 2 ). 

Careful implementations of such algorithms provide acceptable accuracy in computing writhe 
for polygonal curves. But reliably computing the writhe of smooth curves requires another step: 
we must be able to bound the error introduced in approximating a smooth curve by an inscribed 
polygonal curve. The purpose of this paper is to prove: 

X : 

Theorem 1. Suppose C(t) is a simple, closed curve of class C A . We assume C(t) is parametrized 
so that \C'(t)\ > 1, and that we have upper bounds £>i, . . . , on \C'(t)\, . . . \C^ A \t)\. Let C n (t) 
be any n-edge polygonal curve inscribed in C with maximum edge length x and 1/x > 5B 2 . 
If the ribbon formed by joining C n {t) to C{t)for every t is embedded, 

(1) | Wr(C7) - Wr(C n )| <anx 3 + nO(x 4 ). 

where a is a numerical constant less than i?2(5-B| + -B3). 

That is, if the lengths of the edges of C n are approximately constant, the error is bounded by a 
multiple of 1/n 2 . 

The proof is based on Fuller's AWr formula, which gives the difference in writhing number 
between two curves as the spherical area of the ribbon bounded by the curves on S 2 swept out by 
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their unit tangent vectors [|TT|]. (Following Bruce Solomon [|16|], we will refer to such curves as 
tantrices, though they are classically referred to as tangent indicatrices.) 

We begin by defining the writhing number in Section 0. Sections [3| and |] then introduce the 
original form of Fuller's AWr formula. In Sections |5] and |6| we extend Fuller's formula to the case 
where one curve is polygonal and the other is of class C 2 using a natural geometric idea: the tantrix 
of a polygonal curve should be defined to be the chain of geodesic segments on S 2 joining the 
(isolated) tangent vectors of the curve (this was pointed out by Chern in [||]). In the process, we 
discover a surprising fact: the writhe of a polygonal curve is equal to the writhe of any smooth 
curve obtained by carefully rounding off its corners! 

Section [7] contains the remainder of our work: estimating the terms in our improved version of 
the AWr formula to obtain Theorem We test our error bounds in Section [8| by computing the 
writhe of a collection of polygonal curves inscribed in a smooth curve of known writhe. 

The last section contains a discussion of some open problems inspired by the present work. We 
state the most important of them now: Like most of the theory of writhing numbers, the proof of 
our main theorem depends essentially on the fact that C is closed. Can these methods be extended 
to open curves? 



2. Definitions 
The writhing number of a space curve is defined by: 
Definition 2. The writhe of a piecewise differentiable curve C(s) is given by: 
m jwn 1 f C"M x C'(t) ■ (C(s) - C(i)) 

<2) Wr(C) = 5F Lc |C(.)-c(t)P ds di ' 

Definition [Z]is inspired by the Gauss formula for the linking number of two space curves, A(s) 
and B(s) (see Epple for a fascinating discussion of the history of this formula): 

When the two curves A and B become a single curve, their linking number becomes the writhing 
number. This introduces a potential singularity on the diagonal of C x C, but a careful calculation 
shows that the integral still converges. In fact, the integrand of Equation ^ approaches on the 
diagonal of C x C, even when the curve C has a corner. 

From now on, we'll assume that C is simple. With this assumption, another way to look at the 
integral of Definition || is to observe that the integrand is the pullback of the area form on S 2 under 
the Gauss map C x C — > S 2 defined by 

(4) (C(s),C(t))- ' C{s) ~ C[t) 



\C(s)-C(t)\ 



From this perspective, we can see that the (signed) multiplicity of the Gauss map at any point p 
on S 2 is just the number of self-crossings of the projection of C in direction p. 
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3. Fuller's A Wr formula 

Suppose we have a differentiable curve C(t), with unit tangent vector T{t). As we mentioned in 
Section [I], the curve T(t) on the unit sphere is known as the tantrix of C . This curve divides the unit 
sphere into a number of cells. Within each cell, the signed crossing number of the projection of C 
is constant: changing projection directions within the cell amounts to altering the projection of the 
knot by a regular isotopy consisting of Reidemeister moves of type II and III (pictured below). 
Neither of these moves changes the signed crossing number of the knot. 




Figure 1. Changing the projection direction within a cell can only alter the 
diagram by one of these two moves. Neither changes the signed crossing 
number of the diagram, as we can see by counting the + and — markers at 
the crossings of C. 

This observation motivates the idea that the writhe of a closed space curve is related to the 
fraction of the sphere's area enclosed by its tantrix. In 1978, Brock Fuller stated the following: 

Theorem 3. (Fuller's Spherical Area Formula) For any closed space curve C (s) of class C 3 , let A 
be the spherical area enclosed by the tantrix of C. Then 

(5) 1 + Wr(C) = — mod 2. 

Fuller used this formula to conclude that the difference in writhe between two curves X and X 1 
whose tantrices T and T\ are sufficiently close is given by a certain formula, which represents the 
spherical area of the ribbon between T and T\ . 

To be more specific, suppose that X and X\ are simple closed space curves of class C 2 , with 
regular parametrization (that is, parametrized so that X' G and X[ never vanish), and unit tangent 
vectors T and T\. Let F : S 1 x [0, 1] — »■ R 3 be a continuous deformation of X into X±, where 
F(t, A) = X\(t) and the X\ are simple curves of class C 1 , with unit tangent vectors T\(t) contin- 
uous in (t, A). 

Theorem 4. (Fuller's A Wr Formula) IfTiit) and T\(t) are not antipodal for all (t, A), then 

(6) WrpTO - Wr(Xo) = 1- [ ^^/fl, ■ m) + n(t)) dt 

In J c 1 + l (t) ■ li(t) 

We observe that this formula does not require an arc-length parametrization of X and X\. 
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4. Justifying Fuller's interpretion of the AWr formula 

While Fuller stated both these theorems in 1978, he did not provide complete proofs for either. 
The first rigorous proofs of Theorems |3| and |] were given by Aldinger, Tabor, and Klapper [Q] 
in 1995. While these authors proved both theorems as stated, they did not show that the formula 
in Theorem [f] represents the spherical area of the ribbon between T and T\ (in [Q], the right-hand 
side of Equation |^ describes the difference between the twist of two frames on X and X\.) 

In the spirit of their paper, we now justify Fuller's original intuition about Equation ^. 

Proposition 5. Given two curves T (t),Ti(t) : [0, 1] — > S 2 where T (t) and Ti(t) are never an- 
tipodal, the area of the spherical region R bounded by T , T\ and the great circle arcs joining their 
endpoints is given by 

Proof. We let 

u(6,t) = cos6 T (t) + sin0 T x {t), 
and parametrize the region R by 

where 9 ranges from to ir/2. Plugging this parametrization into the area form on S 2 , and using 
the properties of the triple product, we find 

I / Ou du 
dArea = — — x — • u ) d0 A dt. 

\u\ 6 \ou os 

Using the definition of u(9, t), this simplifies to 



dArea = T x T\ 




cos 6 



2 cos 9 sin 9 T • T\ 



sin 9 



2cos9sui9 T • T x )s 



-T[ c\9 A dt. 



Using the formula sin 20 = 2 cos 9 sin 9, and the fact that the definite integrals of each of the 
trigonometric expressions above from to n/2 are equal, we have 



Area(.R) 



T x T x 



it/2 



cos 9 



sin2^T -T 1 ) 3 / 2 



d9 



(Tq + T{) dt. 



This can be solved by the general integration formula 



(8) 



cos 9 



'1 + a sin 20) 3 /2 



d# 



-acos^ — sin 9 



l)Vl + asin20 



which yields the formula in the statement of the Proposition. 



□ 
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5. Extending Fuller's formula to polygonal curves: I 

To measure the difference in writhe between a smooth curve and a polygonal curve inscribed 
in the smooth curve, we must extend Theorem ^ to polygonal curves. To do so, we intend to 
approximate each polygonal curve with a family of smooth curves so that the writhe of the smooth 
curves converges to the writhe of the polygonal curve. 

Examining Definition |2|, it might seem that this result follows from general principles. For 
instance, one might conjecture that Wr was continuous in the C 1 norm on curves, and hope to 
obtain an approximating family using standard techniques. Unfortunately, the situation is not so 
simple; as the example in Figure 2 shows, writhe is not continuous in any C k norm on curves. 
Thus, our proof depends explicitly on the hypothesis that the limit curve is polygonal; it cannot be 
easily extended to the case where the limit curve is merely piecewise C 2 . 




Figure 2. The family of almost-planar curves on the left converge in any C k 
norm to the planar figure eight curve on the right. However, the writhe of 
the curves on the left approaches one, while the writhe of the planar figure 
eight is zero. This shows that writhe is not continuous in any C k norm on 
curves. 

To prepare for the proof, we establish some notation for polygonal curves. Let C(t) be a polyg- 
onal curve with corners at cyclically ordered parameter values t < ti < ■ ■ ■ < t n = t . We let 
T(t) denote the unit tangent to C, and set up the convention that T(£j) will be the tangent vector 

leaving C(ti). 

We now construct a family of smooth curves approximating our polygonal curve. 

Proposition 6. Given an embedded polygonal curve C with corners at to, ... , t n -i,t n = to, there 
exists a family of smooth curves Ci converging pointwise to C with 

1. Ci = C outside a neighborhood of each corner point C '(tj) of radius 

2. Near each corner, the tangent vectors ofCi interpolate between T(t,-_i) and T(tj). 

3. Wr(Cj) -> Wr(C). 

Proof. It is easy to construct a family of Ci — > C obeying conditions (1) and (2) by rounding off 
each corner of C. We claim that this can be done in such a way that the writhe integrand has a 
uniform upper bound on all the Cj. Since condition (1) implies that the C, — > C pointwise in the C 1 
norm, the bounded convergence theorem [[15], p.81] will then yield condition (3). 

Since any pair of adjacent edges is planar, we can choose the Cj so that the region of each Cj 
approximating a pair of adjacent edges is also planar. This means that for some universal e, the 
writhe integrand of each Cj vanishes in an e-neighborhood of the diagonal of Ci x Cj. 

Since C has no self-intersections and the angle at each corner of Cj is positive, the distance 
between any pair of non-adjacent edges of C is bounded below by some constant. Since the Cj 
converge to C pointwise, we may assume the same for the portions of the Cj approximating any 
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pair of disjoint edges. Throwing away finitely many of the Cj if necessary, this means that for any 
5 > 0, there exists a universal lower bound (depending on 5) on the distance between any pair of 
points in Ci x d outside an (^-neighborhood of the diagonal. 

But for any pair of points on Cj, the writhe integrand is bounded above by the inverse square of 
the distance between them. Thus, our lower bound on self-distances yields a universal upper bound 
on the writhe integrand for C and all the Cj outside a 5-neighborhood of the diagonal. Choosing 
5 < e, this completes the proof of the proposition. □ 

6. Extending Fuller's formula to polygonal curves: II 

We now state our extension of Fuller's theorem. Our formula will apply to the following situa- 
tion (c.f. Section |3|): Suppose that X and X\ are simple closed space curves, with X of class C 2 
and X 1 polygonal, with regular parametrization (that is, parametrized so that X' and X[ never 
vanish where they are defined), and unit tangent vectors T and 2\. 

Let F : S 1 x [0, 1] ->• R 3 be a C° deformation of X into X x , where F(t, A) = X x (t), so that 
the X\ are simple curves of class C 1 for A £ [0, 1), with unit tangent vectors T x (t) continuous 
in (t, A). As above, we take the corners of X\ to be at parameter values t , h, . . . ,t n = t . We let 
Ti denote the unit tangent vector to Xi, and let Ti(ij) be the tangent vector leaving Xi(tj). 

Theorem 7. If each corner angle of X\ is strictly greater than n/2, and each Ti(t) and T\(t) are 
at an angle less than it/ 2, then 



Wr(Xi) - Wr(X ) = — V Area #(7^), T (t i+1 ), Tifc)) + Area ATo^T^)^), 

2tx a — ' 

i=l 

where R(T Q (ti) , T (t i+1 ) , Ti (ti)) is the spherical region bounded by geodesies from Ti (ij) to T (ti) 
and T (t i+ i) and the portion of T between ti and t i+ i, AT (t i )T 1 (tj_ 1 )T 1 (t i ) is the spherical 
triangle with these three vertices, and Area represents oriented area on S 2 . 

Proof. Construct a sequence of smooth curves Cj — > Xi using Proposition ^|. For large enough j, 
each of these curves can be homotoped to Xi through a family of simple C 1 curves with a contin- 
uous family of tangent vectors, as in the setup for the statement of this theorem above. 

Joining these homotopies to the homotopy from X\ to X assumed by our hypotheses generates 
a family of (non-smooth) homotopies from the X to each of the Cj. We wish to smooth each 
of these to obtain homotopies from X to Cj which obey the conditions of Fuller's AWr formula 
(Theorem |]). 

We first prove that the tangent vectors of each of the intermediate curves in each homotopy 
from X to Cj are never antipodal to the corresponding tangent vectors Tj of Cj. By hypothesis, 
for each t and A, Z.T\(t), T\(t) < tt/2. On the other hand, since the difference between the tangent 
vectors to X\ at any corner is less than n/2, for large enough j, ZTi(£), Tj(t) < n/2. Putting these 
equations together, we see that ZT\(t), Tj(t) < re, and so these vectors are never antipodal. 

It is easy to smooth the combined homotopy from X to Cj so that each of the intermediate 
curves is of class C 1 while preserving this condition. Since the smoothed homotopy satisfies the 
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hypotheses of Fuller's AWr formula (Theorem^), Proposition |5] tells us that the difference between 
Wr(X ) and Wr(Cj) is the spherical area of the ribbon joining T and Tj. 

For each i, the contribution to the spherical area from the straight part of Cj between U and U + i 
comes from the ribbon between Xi(ij) and the portion of T with t e (ij + 1/j, — As 
j — > oo, this area converges to the area of the ribbon between the portion of T with t 6 (tj, t i+ {) 
and Tiiti). This is the first term in our sum above. 

At each vertex tj of X\, the contribution to our spherical area from the curved part of Cj comes 
from the ribbon between the great circle arc connecting Ti(tj_i) and Ti(t) and a portion of T of 
parameter length 2/j. As j — > oo, the area of this ribbon converges to the area of the spherical 
triangle with vertices T (U), Ti(U), Ti(^_i). This is the second term in our sum above. Figure 3 
shows both these terms on the unit sphere. 



/ / ' \\^7\A T » 



Figure 3. This figure shows the two types of regions in the sum in the 
statement of Theorem [7]. The top (dotted) curve shows the great circle arcs 
joining the tangent vectors T(ti) of the polygonal curve Xi. The bottom 
curve shows the continuous curve of unit tangents T to the smooth curve 
X . The light gray regions show the first terms in the sum, while the dark 
gray spherical triangles show the second terms. 

We have shown that the right-hand side of the statement of the Theorem is equal to the limit 

lim i ^ 00 (Wr(C J ) - Wr(X )). However, by Proposition §, lim^ Wr^) = Wr(Xx). Thus 

(9) lim Wr(C7 i ) - Wr(X ) = Wr(X x ) - Wr(X ), 

which is the left-hand side in the statement of the Theorem. This completes the proof. □ 

We now make a surprising observation: Since the tantrices of the Cj differ as curves on S 2 only 
in parametrization, the area between each of these curves and the tantrix of X is constant. Thus, 
by Fuller's formula, each Cj has the same writhe! And since (by Proposition |6|) these writhing 
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numbers converge to the writhe of X x , each Wr((7j) is equal to Wr(Xx) as well! So we have the 
following corollary: 

Corollary 8. If C n is a polygonal curve, and C is a smooth curve obtained by rounding off the 
corners ofC n under the conditions of Proposition^ then 

(10) Wr(C n ) = Wr(C). 

7. Bounding the AWr Formula 

We now prove our main theorem by finding asymptotic bounds for Fuller's AWr formula. Our 
theorem deals with the following situation: Assume that C(t) is a simple closed curve of class C 4 , 
parametrized so that |C"(£)| > 1. (Given any initial parametrization, this can be accomplished 
by rescaling.) Further, assume we have upper bounds B x , . . . , B A on the norms of the first four 
derivatives of C. 

Let C n (t) be any n-edge polygonal curve inscribed in C. We assume that the maximum edge 
length of C is bounded by x. 

Theorem 9. If the ribbon formed by joining C n (t) to C(t) for every t is embedded, and 1/x > 
5B 2 , 

(11) | Wr(C) - Wr(C n )| <am 3 + nO(x 4 ). 
where a is a numerical constant less than B 2 (5B 2 + B 3 ). 

We make a few comments on this theorem before diving into the proof. First, we observe that if 
the lengths of the edges of C n are all of the same order of magnitude, the difference between the 
writhe of C and the writhe of C n is of order 1/n 2 . 

Next, we observe that the form of our theorem was chosen to be of maximal use in applications. 
In particular, we did not require that C be parametrized by arclength and state our bounds in terms 
of curvature and torsion because in practice it is very difficult to obtain an arc-length parametriza- 
tion of a given curve, while it is comparatively easy to obtain values for the derivative bounds given 
above. 

Last, we discuss the role of the additional hypotheses in the statement above; that the ribbon 
between C and C n be embedded and that 1/x be greater than 5B 2 . Both are intended to exert 
enough control over the approximation to guarantee the existence of a homotopy from C to C n 
obeying the requirements of Theorem [71. 

We can guarantee that C n satisfies the first hypothesis by proving that C n lies in an embedded 
tubular neighborhood of C. Since C is of class C 4 , and has no self-intersections, such a neighbor- 
hood is guaranteed to exist: for a discussion of how to compute the radius of this tube (which is 
known as the thickness of C), see the literature on ropelength of knots (e.g. [ ]12| , [7J, [l3j]). 

Proof. We begin by reparametrizing our curve by arclength. This forces us to recompute our 
bounds for the derivatives of C(t) (a standard computation), arriving at 

(12) |C"0)| = 1, |C"(»| <K:=2B 2 , \C"'(s)\ < T := 2B 3 + 10B 2 2 , 



ON COMPARING THE WRITHE OF A SMOOTH CURVE TO THE WRITHE OF AN INSCRIBED POLYGON 



9 



while C^'(s) is again bounded above. To remind ourselves of the connection between these 
bounds and the curvature and torsion of our curve, we will refer to the bound for the second 
derivative as K, and the bound for the third derivative as T. Further, we note that the curvature 
k(s) of our curve is bounded above by K, and that our hypotheses imply that l/x > (5/2)K. 

We also establish the convention that the corners of C n are at parameter values cyclically ordered 
as so, . . . , s n _i, s n = sq. 

By smoothing the linear interpolation between C and C n , we can construct a homotopy between 
C and C n according to the conditions of Theorem [7| as long as: 

1. the ribbon joining C to C n is embedded, 

2. the angle at each corner of C n is at least tt/2, 

3. the angle between T(s) and T n (s) is at most %/2 for any s. 

Borrowing from Lemma [12] (below), we see that our assumption that l/x > (5/2) Hf is enough 
to bound the angle in (3) by 0.20402 < 7r/4. At any corner Sj, the same Lemma implies that the 
corner angle is the supplement of at most twice 0.20402, so this is enough to ensure that condition 
(2) holds as well. 

Theorem [7] now tells us that 

1 n 

(13) |Wr(C)-Wr(C„)| < — V |Areai? (T (s t ) , T (s i+1 ) , T n ( Si ))| 

ZTT * — ' 

1=1 

+ | Area AT (sj) T n (si_i)T„ | , 

where the first term is the area of the spherical region bounded by the geodesies from T n (si) 
to T(si) and T(sj + i) and the portion of T between Sj and s i+ i, and the second term is the area 
of the spherical triangle. Our job now is to estimate the areas of these regions. To do so, we first 
invoke Taylor's Theorem, in the form commonly used in numerical analysis (c.f. [|3|], Thm.1.4): 

Theorem 10. (Taylor's Theorem) Suppose C(s) is a curve of class C 4 , with fourth derivative 
bounded by B' 4 . Then (choosing coordinates so that C(0) is at the origin), 

(14) C(s) = sC"(0) + £ C"(0) + ^ C"'(0) + R,(s), 

2 o 

where \R A (s)\ < s 4 ^. 

We will use this expression for C(s) frequently in our work below. 
Lemma 11. For any s, we have 

(15) | S _|C( S )||<^| S 3 | + T L| S 5 |, and |C(s)|<|4 
Further, for any edge ofC n , the difference \s i+ i — Sj| is at most 1.01 x. 

Proof. We assume without loss of generality that s is positive. By Schur's lemma ([§]), since the 
curvature of C is bounded above by K, \C(x)\ is at least the length of a chord across an arc of 
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length s on a circle of radius l/K, or (2/K) sm(K/2)s. This means that we have 

2 K K 2 3 r, / N 

— sinys = s-— s +i? 5 (s), 

where R$(s) is the term of order s 5 which comes from the usual Taylor expansion of sins. In 
particular, 

2 K 
\s-\C{s)\\ < s-— sin—s 



K 2 



24 

where Rq(s) < -^s 5 . The upper bound on |C(s)| comes from the fact that C is unit-speed. 

The second statement is another Schur's lemma calculation; this time invoking our hypothesis 
that x > (h/2)K and observing that 1.01 sin?/ > y for y between and 1/5. □ 

We will also need an upper bound on the angle between T(s) and C n (s). 

Lemma 12. The angle between the tangent vector T(s) and the corresponding tangent vector 
T n (s) to C n is bounded above by 

(16) ZT(s)T n (s) < 0.510057^. 

Proof. Assume that s is between Sj and Si + %. Then 



But we have 



and for any t, we have 



C(s l+1 )-C( Sl )= I l+1 T(t)dt, 



Si 



T(t) = T{s) + J T\u) du. 



This means that 

(18) [C(s i+1 ) - C{ Si )\ x T{s) = [* +1 T(t)xT(s)dt 

J Si 

(19) = f 1 f T\u) xT{s)dudt. 

J S{ J s 

Since \T'(u) x T(s)\ < \T'(u)\\T(s)\ < k(u) < K, and s is between Sj and s i+1 , a small compu- 
tation reveals that this integral is bounded by v( s i+i — s «) 2 - 

Since the length \C(s i+ i) — C(sj)| is bounded below by (1/1.01) — Sj) by Lemma [TT], we 
get 

(20) sin ZT{t)T n {t) < ~^-Kx. 
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Since 1/x > (5/2)K, this is always bounded above by 1.01/5, and so 

1 01 2 

(21) ZT(t)T n (t) < -^-Kx. 



□ 



We are now ready to embark on the main work of the proof: estimating the areas in Equation [T3|. 
We begin with the first term: the area bounded by the portion of T(s) between Sj and s i+1 , together 
with the great circle arcs joining T(s$) and T(s i+1 ) to T n (sj). Without loss of generality, we may 
assume that i = 0, that s = 0, and that C(0) = 0, and apply the Taylor expansion of Equation [TJ] 
to C at 0. Our strategy is to prove that this region is contained in a neighborhood of the great circle 
arc joining T(0) and T(s 1 ). Suppose s is between and s\. We want to bound the height of T(s) 
above the T(0), T(s{) plane, or 

(22) k{s) " 1^x^)1 • 

First, we have 

C'( Sl ) = C'(0)+s 1 C"(0) + ^C'"(0)+R 3 (s 1 ). 

C'(s) = C'(0) + sC"(0) + jC"'(0) + R 3 (s). 

Using the triple product identities, we can rewrite h(s) in terms of the inner product of C"(0) and 
the cross product of these vectors. For the triple product, we get 



(23) 



sjs s 2 s 1 



C"(0) • C"(0) x C'"(0) + C"(0) • [R 3 ( Sl ) x C'{s) + C'(s x ) x R 3 (s)} . 



Expanding the last term, we see that is the sum of a term of order Sis 3 and a term of order ssf . 
Thus, to leading order, the norm of the entire triple product is bounded above by 

(24) ^ <H --= 2\CmxC> M \ S ° + °^ 

since s G [0, si]. We now consider the height of T n (0) above the T(0), T(s\) plane. Since T n (Q) 
is the normalization of C(s\) — (7(0) = C(si), this height is given by 

C( Sl ) C'(0) x C"( Sl ) 



|C( Sl )| |C7'(0)xC7'( Sl )| 



(25) 

As before, we get 

(26) t7'(0) x C\ Sl ) = Sl C'(0) x C7"(0) + ^ C7'(0) x C""(0) + 0(s?). 

Taking the dot product with the Taylor expansion of C(si), we get only terms of order 0(sf) and 
higher. Thus, to leading order, this region is contained in a rectangle based on the great circle arc 
joining C"(0) and C'(si) of height H. We now estimate the area of this rectangle. 
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First, we note that the length of the great circle joining C'(0) and C'(si) is given by the angle 9 
between C"(0) and C'(si). Since si < 1.01a; by Lemma [TTJ, this length is bounded above by 
1.01 Kx, which is less than 0.404 by our hypotheses on x. Since H is small compared to s, we 
may assume that the entire rectangle is contained within a spherical disk of radius 0.5. 

We project the rectangle to the plane by central projection: this map is increasing on lengths and 
areas, and increases length by at most a factor of 1.01. The area of the rectangle in the plane is 
overestimated by the product 1.01 9H. On the other hand, we have \C'(0) x C'(si) \ = sin 9. And 
for 9 < 0.404, 1.02 sin 9 > 9. Keeping track of the various constants involved, and using the fact 
that si < 1.01 x again, the area of this spherical region is overestimated by 

(27) Area R(T( Si ),T{s i+1 ),T n ( Si )) < KTx 3 + 0(x 4 ), 

We now turn to the second term in the Equation [Of the area of the spherical triangle bounded 
by T(sj), X„(sj_i) and T n (si). Without loss of generality we assume that i — 1, that si = 0, and 
that C(0) = 0, and we expand C around using Equation |T3]. We wish to compute 



(28) Ah. ^(-^4,^4,^(0) 
,|C(s )| \C(s 2 )\ 



If we factor out \/\C(sq) ||C(s 2 )|, we are left with the norm of the cross product of two terms: 

C(s )-\C(s )\C'(0) = (s o -\C(s o )\)C'(0) + S fC"(0) + ^C"'(0) + R A (s o ) 
C(s 2 )-\C(s 2 )\C'(0) = ( S2 -|C( S2 )|)C"(0) + |c"(0) + |c w (0) + J R 4 ( S2 ). 

2 D 



Using Lemma |TT], we see that \s — \C(s)\ \ < (K 2 /24)s 3 + 0(s 5 ), and we see that the leading term 
of this expression contains fifth powers of of s and s 2 , and is bounded by: 

K 3 KT 
48 + ~12 



(29) 44 [ ~— + "-T7T I (SO + S 2 



However, we must still divide by |C(s ) | |C(s 2 ) |- By Lemma [TTi we see that the ratios s /\C(s )\ 
and s 2 /|C(s2)| are bounded above by 1.01. Thus, using the same Lemma to conclude that s 2 and s 
are less than 1.01 x, and making a central projection argument as before, we are left with 

(30) AreaA^^),^^)^^)) < ^ x 3 + 0(x 4 ). 

Summing over i, and dividing by 2n, then writing K and T in terms of B 2 and B 3 , we obtain the 
statement of the theorem. Note that we have overestimated the numerical constants to simplify the 
resulting formula. □ 

If a curve has a small region of high curvature, and larger regions of low curvature, it may be 
desirable to approximate the curve more carefully in the regions of high curvature in order to save 
time in the computation of writhe. Since our error bound is additive along the curve, these methods 
are well suited to this case. We have 
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Corollary 13. Suppose C is a C 4 curve and C n is a curve inscribed in C so that C and C n obey 
the hypotheses of Theorem ^. 

If C and C n are divided into regions Ri, each containing rii edges which are bounded above in 
length by Xi, and so that the derivatives of C are bounded by Bu, . . . , B& and l/x{ > 5B2i, then 

\Wv(C)-Wv(C p )\<J2^ n ^i+ n i°( x t)- 

i 

where each on is a numerical constant less than B-iii^B^ + B^). 

We make one more observation: 

Proposition 14. Let C be a simple, closed space curve of class C 2 , and C p be a polygonal approx- 
imating curve as in Theorem ^ or Corollary [D[ 

If the arc joining the endpoints of a sequence of n edges of C p is planar, then the n — 2 edges 
interior to this region contribute nothing to the error bound in the Theorem. 

In particular, this means that the derivative bounds in both statements can be taken to be bounds 
on the derivatives of the non-planar regions of the curve C. 

Proof. On these edges, the tantrix of the smooth curve and the polygonal curve parametrize the 
same great circle arc on S 2 . Thus, the ribbon between these curves has zero area. □ 

8. Example Computations 

We are now prepared to test Theorem [j] by computing the writhing numbers of various polyg- 
onal approximations of a smooth curve, and comparing the results to the exact writhe of the 
smooth curve. To control the numerical error introduced in these calculations, all of these compu- 
tations were performed using an arbitrary-precision implementation of Banchoff 's formula for the 
writhing number of a polygonal curve. The initial runs were performed with 45 decimal digits of 
precision. They were checked against runs performed with 54 digits of precision. Since the results 
agreed, we feel confident that roundoff error does not affect the computations reported on below. 

The curve whose writhe we computed is an example of Fuller[JTT]]: 




Figure 4. This example of Fuller's "closed helix" is composed of 3 turns of 
a helix of radius 1 with pitch angle 0.33, with ends joined by a planar curve. 

Using Theorem |3], and the Calugareanu-White formula, it is easy to see that the writhe of this 
curve is 3(1 — sin 0.33) ~ 2.0278709. After all, the area enclosed by the tantrix of this curve C is 
that of a hemisphere, plus 3 enclosures of a spherical cap of radius n/2 — 0.33. Thus the writhe 
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of the curve is equal to 1 — sin 0.33 mod 2. To complete the computation, one sets up a frame on 
the curve, and computes its twist and linking number. (Details for this computation can be found 
in UIJ.) 

We now take a series of polygonal approximations to C, and compare the difference between 
their writhing numbers and the writhe of C to the bounds of Theorem ^. 

We begin by finding bounds on the derivatives of C and the edge length of our approximations. 
By Proposition [13], it suffices to find derivative bounds for the helical region of C . Since the 
helix has unit radius, both B 2 and B 3 can be taken to be one. The curve is parametrized so that 
\C'(s)\ > I. 

Here are the results of computing writhe with various numbers of edges: 



n 


Wr(C n ) 


|Wr(C n )-Wr(C)| 


X 


anx 3 


100 


2.00541 


0.02246 


0.506 


77.73 


250 


2.02434 


0.00353 


0.203 


12.55 


500 


2.02697 


0.0009 


0.101 


3.09 


1000 


2.02763 


0.00024 


0.051 


0.786 



It is worth examing a graph of these results. 

log anx .ie3 

Ae2 
I. 

.1 

log|Wr(C„)-Wr((7)| ^ 

.le-2 



le3 .2e3 .4e3 .7e3 .le4 




Figure 5. This graph shows a log-log plot of the actual error in comput- 
ing the writhing number for one of Fuller's "closed helices" with various 
numbers of edges (lower solid line), together with our error bounds (upper 
dotted line). The fact that the lines are parallel shows that the convergence 
is of order n 2 , as predicted by Theorem |9|. 



9. Further Directions 

In this paper, we have given a set of asymptotic error bounds which allow us to compute the 
writhe of a closed space curve with defined accuracy by computing the writhe of a polygonal 
approximation to this curve. The example we computed in Section || shows that our bounds are of 
the right order of magnitude: roughly speaking, the writhe converges quadratically in the number 
of edges of the approximation. Our work leaves open several directions for further inquiry. 
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First, it is puzzling that our approximation theorem for curves with corners (Proposition |6|) 
should depend on the hypothesis that the limit curve is polygonal. To sketch an extension of this 
result, we recall a definition from Chern ([§]]): 

Definition 15. The tantrix of a piecewise C 1 curve C(s) with positive corner angles is the image 
ofT(s) on the unit sphere, together with the great circle arcs joining the pairs of tangent vectors 
at each corner of the curve. 

We note that our Theorem [7] shows that Fuller's A Wr formula holds for polygonal curves with 
the definition of tantrix extended as above. We further suspect that: 

Conjecture 16. Fuller's Spherical Area formula (Theorem [|j and Fuller's AWr formula (Theo- 
rem ^) hold for piecewise C 2 curves with the extended definition of tantrix given by Definition [73|. 

The proofs of both of these theorems depend on the Calugareanu-White formula, which only 
applies to closed curves. Thus all of our results are restricted to closed curves. This leaves open a 
much more important problem: 

Problem 17. Extend all these theorems (the Calugareanu-White formula, Fuller's spherical area 
formula, and Fuller's AWr formula) to open curves. 

In particular, extending the results of this paper to open curves would be useful for applications 
in biology, where the curves of interest are not neccesarily closed. We note that while Fuller's 
AWr formula makes sense for open curves, computational examples show that it does not give the 
correct answer: boundary terms must be added to account for the ends of the curves. 

10. Acknowledgements 

I am grateful to Herbert Edelsbrunner, Herman Gluck, Issac Klapper, John Maddocks, Kath- 
leen Rogers, and David Swigon, among many others, for fruitful conversations. This work was 
supported by both a National Defense Science and Engineering Graduate Fellowship and an NSF 
Postdoctoral Research Fellowship (DMS-99-02397). 

REFERENCES 

[1] Pankaj K. Agarwal, Herbert Edelsbrunner, and Yusu Wang. Computing the writhing number of a polygonal knot. 

In Proc. 13th Ann. ACM-SIAM Sympos. Discrete Alg., 2002. 
[2] J. Aldinger, I. Klapper, and M. Tabor. Formulae for the calculation and estimation of writhe. J. Knot Theory and 

Its Ramifications, 4(3):343-372, 1995. 
[3] Kendall E. Atkinson. An introduction to numerical analysis. John Wiley & Sons Inc., New York, second edition, 

1989. 

[4] T. Banchoff. Self linking numbers of space polygons. Indiana Univ. Math. J., 25(12): 1 171-1 188, 1976. 
[5] George Calugareanu. L'integrale de Gauss et l'analyse des nceuds tridimensionnels. Rev. Math. Pures Appi, 
4:5-20, 1959. 

[6] George Calugareanu. Sur les classes d'isotopie des nceuds tridimensionnels et leurs invariants. Czechoslovak 

Math. J., ll(86):588-625, 1961. 
[7] Jason Cantarella, Robert B. Kusner, and John M. Sullivan. On the minimum ropelength of knots and links. 

Preprint. arXiv:math.GT/0103224. 



16 



JASON CANTARELLA 



[8] S. S. Chern. Curves and surfaces in euclidean space. In S. S. Chern, editor, Studies in Global Geometry and 

Analysis, pages 16-56. Math. Assoc. Amer., 1967. 
[9] D. Cimasoni. Computing the writhe of a knot. J. Knot Theory Ramifications, 10:387-295, 2001. 
[10] Moritz Epple. Orbits of asteroids, a braid, and the first link invariant. Math. Intelligencer, 20(l):45-52, 1998. 
[11] F. Brock Fuller. Decomposition of the linking number of a closed ribbon: A problem from molecular biology. 

Proc. Nat. Acad. Sci. (USA), 75(8):3557-3561, 1978. 
[12] Oscar Gonzalez and John H. Maddocks. Global curvature, thickness, and the ideal shapes of knots. Proc. Nat. 

Acad. Sci. (USA), 96:4769-4773, 1999. 
[13] Richard A. Litherland, Jon Simon, Oguz Durumeric, and Eric Rawdon. Thickness of knots. Topol. Appl, 

91(3):233-244, 1999. 

[14] William F. Pohl. DNA and differential geometry. Math. Intelligencer, 3(l):20-27, 1980/81. 
[15] Halsey L. Royden. Real analysis. Macmillan Publishing Company, New York, third edition, 1988. 
[16] Bruce Solomon. Tantrices of spherical curves. Amer. Math. Monthly, 103(l):30-39, 1996. 
[17] De Witt Sumners. Lifting the curtain: using topology to probe the hidden action of enzymes. Match, 34:51-76, 
1996. 

[18] James H. White. Self-linking and the Gauss integral in higher dimensions. Amer. J. Math., 91:693-728, 1969. 



